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Abstract 

The isospectrality problem is studied for the operator of the spherical hydromagnetic a 2 — dynamo. It is 
shown that this operator is formally pseudo-Hermitian (J— symmetric) and lives in a Krein space. Based 
on the J— symmetry, an operator intertwining Ansatz with first-order differential intertwining operators is 
tested for its compatibility with the structure of the a 2 — dynamo operator matrix. An intrinsic structural 
inconsistency is obtained in the set of associated matrix Riccati equations. This inconsistency is interpreted 
as a no-go theorem which forbids the construction of isospectral a 2 — dynamo operator classes with the help 
of first-order differential intertwining operators. 



1 Introduction 

The magnetic fields of stars and planets are generated by the homogeneous dynamo effect in moving electrically 
conducting fluids 1 . This effect is explained within the framework of magneto-hydrodynamics (MHD), but 
its experimental demonstration was missing until recently. In 1999, the first successful dynamo experiments 
in Riga and Karlsruhe 2 opened up a new way for the laboratory investigation of homogeneous dynamos. 
In connection with the data analysis for the existing experiments and the design of new dynamo experiments 
there is a growing interest in the spectral properties of dynamos. Of particular interest is the question whether 
isospectral dynamos can exist. First numerical results on this topic were published in [13 E], but rigorous results 
are still missing. 

As a step towards clarification of this issue, we study in the present paper the question whether operator 
intertwining techniques from quantum mechanics (QM) can be adopted to MHD dynamo models. In case of 
an affirmative answer we would obtain an efficient tool for constructing isospectral classes of MHD dynamo 
operators. Otherwise we would get a no-go theorem which would forbid a straight analogy with quantum 
mechanical models. 

Let us start by recalling some essentials of operator intertwining transformations in QM 5 . Two operators 
Hq and H\ are said to be intertwined if there exist operators A + and A- so that 

H X A + = A + H , A-Hx = HqA_. (1) 

For the corresponding eigenfunctions <j)Q and cf>i holds, up to normalization, 

(f>i = A + 4>q , A-4>i = 4>q 

and the operators H and Hi are isospectral, except for those states that are annihilated by A + or A_ . In 
the case of one-dimensional Schrodinger operators Hq = p 2 + Vo(x) and Hi = p 2 + Vi(x) with the momentum 
operator given as p = —id x , the intertwining operators can be chosen as first order differential operators 

A+=A:=tp + f, A_=Ai = -ip+f. (2) 

Structural compatibility of Hq and Hi with the intertwining relations requires that the function f(x) and 
the potentials Vq (x) , Vi (x) are connected by the consistency conditions 

Vi = V Q + 2f 
-f + f = Vq-E (3) 
f' + f = Vi-E, (4) 
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where the prime denotes differentiation with respect to x; and E is a constant of integration. Linearization 
of the Riccati differential equations |J5J, I0J shows that this constant can be interpreted as eigenvalue of the 
Schrodinger operators Hq and Hi 

HoXo=E Xo for / = -^, (5) 

Xo 

H lX i=E Xl for f = ^-, (6) 

Xi 

where xo an d Xi are formal, and not necessarily normalized eigenfunctions of Hq and Hi, respectively. They 
are connected by the product relation 

XoXi = c (7) 

with c a non- vanishing constant. It is straightforward to verify that the shifted Schrodinger operators are 
factorizable in terms of the intertwining operators 

H -E = A*A, H 1 -E = AAK 

First-order differential intertwining transformations of type are known as Darboux transformations [H] and 
are widely used to generate isospectral operator classes from given operators with known spectra |SJ [7| |S] . In 
particular, intertwining constructions are a basic ingredient of super-symmetric quantum mechanical models 
03 and their generalizations to pseudo-supersymmetric systems ^| m ■ As it was demonstrated in ^2] > 
a double-intertwining (double commutation) method can provide a tool for inserting additional eigenvalues in 
spectral gaps of given background Schrodinger and Jacobi operators. 

Motivated by the large number of exact results on isospectral classes obtained by operator intertwining 
constructions, it is natural to investigate whether MHD dynamo operators are also suitable for this technique. 
For this purpose we study in the present paper the simplest mean-field MHD dynamo configuration — the 
spherical a 2 — dynamo 1 . In terms of the radial momentum operator p = —i(d r + 1/r) the 2x2 operator 
matrix of the a 2 — dynamo is given as 1 

-P 2 ~ ^ a(r) 
Hi[a]=\ | (8) 

pa(r)p + a(r) l -^l -p 2 - 



and lives on the domain 

V{Hi[a]) :=U>= ^ ^ J : ^eH=H®H, H = L 2 (Q, r 2 dr), 
ft = [0,1], V(1) = 0, r7Hr)U ^0} 

in the Hubert space Ti. It describes the coupled I— modes of the poloidal and toroidal magnetic field components 
in a mean-field dynamo model with helical turbulence function a(r). The function a(r) does not depend on I 
and we assume that it is real-valued, positive definite, bounded, and sufficiently smooth on fi: Im a = 0, < 
a < ci < oo, a € C 4 (f2). The idealized boundary condition ijj(r = 1) = corresponds to a super-conducting 
spherical boundary shell and is chosen here to ensure simplicity of the subsequent considerations |13j . For more 
realistic models with close relation to stellar dynamos, the spherical fluid configurations confined to r < 1 can 

be assumed as embedded in empty space. The boundary condition should then be replaced by Bitp = 

with Bi = diag[<9 r + (I + l)/r, 1] (see e.g. pQ) what will require a more general approach than that presented 
in the present paper. 

Exploring the fundamental symmetry of the a 2 — dynamo operator matrix we find in section [21 that Hi[a] 
acts as a symmetric operator on the Hilbert space H when this is endowed with an indefinite metric J. I.e., the 
a 2 — dynamo operator matrix is a J— symmetric (formally J— self-adjoint) operator 

Hi[a] = Hf[a] := JH}[a)J 

living in a Krein space K. = Hj |14|. J— self-adjointness is a natural property of operators from different 
fields of physics. Examples are, e.g., the super-symmetric Dirac operator |15| . PT— symmetric non-Hermitian 

X A brief outline of the derivation of the a 2 — dynamo operator matrix from the MHD mean-field induction equation can be found 
in Appendix 1X1 
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Hamiltonians in QM ^3 ^5] as well as the Wheeler-DeWitt operator for a cosmological Friedman-Robertson- 
Walker model coupled to a real massive scalar field JUj. Since the recent paper series ^3 El \U\ of A. 
Mostafazadeh on non-Hermitian operators with real spectra, J— self-adjoint operators are also known as pseudo- 
Hcrmitian operators. 

In analogy with the simple quantum mechanical model described above, we base our isospectrality analysis 
on an intertwining Ansatz for two a 2 — dynamo operators with helical turbulence functions ao(r), otx(r) 

H lo [a ] -EI= -AM, H h [ ai ] - EI = -MA 

and intertwining operator matrices A, A* that are first-order differential operators 

A := iR(r)p + Q(r), i J := -ipR?(r) + Q B (r) . 

This Ansatz leads to a set of six consistency conditions on the matrices R(r) and Q(r) which are studied in 
section 03 It is shown that one pair of conditions fixes the structure of R(r) in terms of the helical turbulence 
functions ao{r) and ai(r). A second pair is equivalent to the symmetry relations B = U = on the 
matrix functions 

B := R*Q 

U := R[Q^ - = RBR- 1 - R(R*)' 

and can be regarded as an implicit consequence of the J— pseudo-Hermiticity of the operator matrices Hi [ao] 

and H^lai]. (The prime denotes the derivative with respect to r.) The remaining two conditions can be 
transformed into a pair of coupled matrix Riccati differential equations (MREs) on B and U . 

The consistency of the six conditions is analyzed in section 0] with the help of a step-by-step reduction 
of their complexity. First, we conclude from the limiting behavior of the MREs for r — > that the angular 
mode numbers Iq and l\ in the two dynamo operator matrices should be connected by the incremental relation 
h = io + 1- Then we use the J— symmetry of B to derive from the coupled MREs a system of coupled non- linear 
ordinary differential equations (ODEs) involving the helical turbulence functions ao(r) and ai(r). Analyzing 
these ODEs we are able to show the existence of an inherent contradiction between them. As implication, we 
arrive at a no-go theorem which states that the six consistency conditions cannot be fulfilled simultaneously 
and that, hence, the structure of the a 2 — dynamo operator matrices is not suitable for an operator intertwining 
technique based on an Ansatz with first-order differential intertwining operators. 

In the concluding section [S] we briefly discuss some other methods which could be useful for studying 
isospectrality issues of the dynamo operator matrix and which possibly could provide a technique to construct 
classes of isospectral spherical a 2 — dynamo operators. 



2 J— Symmetry of the dynamo operator matrix 

In this section, we study the fundamental symmetry J of the a 2 — dynamo operator matrix (jSJ) what allows us 
to choose an appropriate Ansatz for the intertwining operators A + and A-. 
We start our consideration by introducing the auxiliary operator 

Q[a\ := pap + a — — — 

defined on the domain 

V{Q) = {0 : <(> € H = L 2 (n, r 2 dr), 0(1) = 0, r^(r)\ r ^ -> 0} 

in the Hilbert space H. The operator Q[a] is a formally self-adjoint 2 singular differential operator Q = 
which acts as symmetric operator on TL. In terms of Q[a] the dynamo operator matrix and its formal adjoint 

_( -Q[i] " \ ivtui_ f -QW <2M 



H M=\ Q[a] -Q[l] h H iM-\ a - Q[1] 



2 In the subsequent compatibility analysis of the operator intertwining construction we restrict our attention to symmetric 
(formally self-adjoint) operators. For simplicity, we leave questions of self-adjoint extensions and corresponding generalized boundary 
conditions 18 19] for the bi-component functions ip aside. 
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so that the fundamental (canonical) symmetry can be obtained as 











1 1 ) 

1 1 ) 


v = 







H l [a]=Hl[a]:=JH}[a]J, J = ( ° J ). (!)) 

Diagonalizing the matrix J 

t 1 
J — ► X] = S JS, S = —pz 

we see that ///[a] is equivalent to the operator matrix 

tr r l cT/rr in 1 ( Q[a-2]+a — Q[a]+a 

ff,[a] = S H t [a}S = - ^[a] ^ Q g[ _ a L _ J 2 ] - a 

with the property flj[a] = r/_ff; [aj^r/. The fundamental rj— symmetry of the operator matrix fij[a] implies that 
T>(Hi) could be endowed with the indefinite metric r\ so that Hi[a] becomes a symmetric operator on V(Hi). 
Due to the invariance of the signature under the transformation S the domain T>(Hi) can also be endowed with 
a natural indefinite inner product [■, -]j defined by the metric J 

[x,y]j := (x,Jy), x,y e H = H ®H, 

where (•,•) denotes the usual inner (scalar) product in the Hilbert space TL. This means that .Hz [a] is a 
J— symmetric operator which acts as symmetric operator in a Krein space 3 K, = Tij 

H t x,y ^ = x,Hfy 

From its operator-matrix representation ® we see that J is self-adjoint, involutory and unitary 

J f = J, J 2 = /, J" 1 = J 1 

so that Hz [a] is a J— pseudo-Hermitian operator in the sense of [TU1 ITT1 IT7] . 

The eigenvalues of J— pseudo-Hermitian operators are known ^3 ^] |20| to be either real or to come in 
complex-conjugate pairs. Here we illustrate this property by passing from the eigenvalue problem for the linear 
operator pencil 

Li[a, \}ip := (# ( [a]-A)v = 

to the eigenvalue problem Li[a, X]ipi — for the associated quadratic operator pencil Li[a, A]. This pencil can 
be derived explicitly from the Ansatz 



^~ ( i[Q(i) + A]£ ) 



with a(r) ^ 0. As result we obtain 



The operators 



Li{a,X]^i = UQ[l] + X]^[Q[l] + X]-Q[a}\ r { =n 
= (A 2 X 2 + AxX + Ao)ipi = 0. 



A :=Q[1]-Q[l]-Q[a], A x := Q[l}~ + -Q[l], A 2 
a a a 



are formally self-adjoint on T>(Q) so that the functionals a,j :— (Ajipi, ipx), j = 1, 2, 3 are real-valued: Im a,j = 0. 
From the quadratic equation 

(Li[a, X]4>i,ipi) = a 2 X 2 + aiX + a Q = 



3 For surveys on operators in Krein spaces (Hilbert spaces with additional indefinite inner product structures) we refer to the 
mathematical literature 1141 1201 . 



4 



we conclude that the eigenvalues of the J— pseudo-Hermitian dynamo operator matrix Hi and its associated 
pencil Li occur as eigenvalue pairs |21| 

1 



A± = — I -01 ± sj a\ - 4a Q a 2 J . 

Obviously, the sign of the discriminant A :— a\ — 4aoa2 defines whether A± are both real or pairwise complex 
conjugate. The transition from real eigenvalues A± to complex ones occurs at A = where the eigenvalue 
becomes two- fold degenerate A+ = A_ = Ao = — This general behavior of A± confirms the results of 
numerical simulations [2| , which showed that a scaling of the helical turbulence function a leads to a pairwise 
intersection of real eigenvalue branches of Hi and a transition at the intersection points to a pair of complex 
conjugate eigenvalues. 

We note that at the two-fold degenerate points Ao of the spectrum with A = a Jordan-Keldysh chain 22 
exists for the linear operator pencil 

L ; (A )V = 0, Li(\ ) X = i> 
as well as for the quadratic operator pencil 

Li(Ao)Vi = , MA )xi + 9aMWx = 0. 
Both are built up from eigenvectors ip, ipi and associated vectors \i Xii respectively. 

3 Consistency conditions and matrix Riccati equations 

The fundamental J— symmetry ( J— pseudo-Hermiticity) of the a 2 — dynamo operator matrix provides a natural 
Ansatz for an intertwining construction which respects this symmetry: 

H to [a ] - EI = -AM, H h [ax] - EI = -i J A. (10) 

In general, the operator matrix A could be an nth-order differential operator of the form 

n 
k=l 

with 2x2 matrices Rk(r) and Q{r) as coefficients. For simplicity, we restrict our attention in the present paper 
to the first-order differential operator 

A = iR(r)p + Q{r) (11) 
with J— adjoint A$ — —ipR$(r) + Q"(r). Here we define the jj— operation for a given 2x2 matrix C as 

C s = JC^J = JC* T J. 

Asterisk and superscript " T " denote complex conjugation and transposition, respectively. 
Let us introduce the abbreviations 

K i := I — cxqact- 



\/r ts fo,iOo,i + l) . rp T 

A/04 : = Ko,i 5 \-EI-ao,i<Tj 



with the nilpotent matrices <j± defined as a + := f g )' a ~ ' = ( 1 J' sm ft e< i a 2 — dynamo 

operator matrices in (|10|) take then the short form 

H la i [00,1] - EI = -pK 0A p - M ,i. (12) 

Substituting (jl 1|> and l|12|l into the intertwining Ansatz l|10|) . making use of commutation relations like [p, R(r)] = 
—iR'(r) and equating the coefficient matrices of the p 2 , p, I terms we obtain the following six consistency 
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conditions 

H lo : p 2 : RR* = K a (13) 

p: RQ i - QR i - + R'B i = (14) 

/ : QQ S - R{R S )" + i?(Q 8 )' - Q(R S )' = M (15) 

H h : p 2 : R}R = K x (16) 

p: -R % Q + Q»i? = (17) 

I : Q S Q - {R S Q)' = Mi. (18) 

For a successful intertwining construction these matrix equations should be fulfilled simultaneously. So, the 
main task consists in finding explicit solution sets for ((13(1 - 1)18(1. Alternatively, we should obtain intrinsic 
contradictions within this equation system which could be interpreted as a no-go theorem forbidding this 
construction for a 2 — dynamo operator matrices. 

We start our analysis with eqs. {Eft and JEJ. From the tautologies RR}R = RR}R and RpRR? = R}RR} 
follows 

RKi = K R, = R^Kq 



what with 
yields 

Hence, we can set 
Using this and 1(19(1 in 

we find 



R 



rii 



( rn 


Til \ 


1 


R* = 


( rh 




\ Til 


T22 J 






\ T*21 


r*n ) 




= 0, 


Oil 


_ Tu 






Til 


«0 


T22 


r* 






Til = 


\r22 




T21 = 


|r 2 i|e 1 ^ 



(19) 



= K Q = ( 1 9 ) , R*R = Ki=( 1 







where the phases 7 and e are still undefined. 

As next step we analyze Eqs. 1(14(1 and 1(17(1 . It is easily seen that defining the matrices 

U := R[Q S - {R s )'] , B:=R*Q (21) 
these equations are equivalent to the J— symmetry relations 

u = u\ b = bK 

Due to the different symmetry content of B and Q it is natural to consider B as primary structural element of 
the intertwining construction, and Q as a secondary one. So, we perform our subsequent investigation in terms 
of B and R. Explicitly, the J— symmetry is realized by the matrix structure 

B =( h+ b : h h \)> ^ = °> * = l,-.,4. (22) 
Furthermore, we exclude Q from ((21(1 to obtain 

U = RBR- 1 - R(R*)'. (23) 



G 



Introducing the notation TV := R 1 R' and substituting ((23(1 into the symmetry relation U — U* yields the 
additional constraint 

[B,Ki x ] = AT* - N. (24) 

From eq. 1(201) we find 



N = i-y'I + 



f 



1 fa' 



a 



-q 
/ q 



Ol 

' 2 



(1 + itane) + i - 

ao cos 2 



(25) 
(26) 



so that (|24|l transforms to 







2q 
/*-/ -2q 



— 2ib^ —hi 

Finally, we arrive at the following restrictions on the phase 7 and the components 62 and 64 of the matrix B: 



i = 0, 



2q 



Im / 

ai 



1 (oL 



• tane - 



a 



(27) 



Summarizing the implications of the first four consistency conditions we see that they are free of intrinsic 
contradictions. From the initially eight arbitrary complex- valued functions contained in the matrices R and Q 
only the three real- valued functions (61, 63, e) are still undefined. Together with the helical turbulence functions 
(ao, cti) and the constants (7, E, Iq, h) £ K 2 x 1? + we expect them to be highly fine-tuned by the remaining two 
consistency conditions (|15|l and l|18|) . 

Let us study these conditions now. Making use of the definitions of U and B in 1(21(1 . their implications 



q* - (B}y 
(q»)' - {R*y 

Q 



-R- l R'R- l U + R- X U' 
(R^B 



(28) 



(29) 



and setting at the end RR^ = Ao, R^R = K\ according to eqs. \Y,\\. we find that (|15|l and 118(1 transform 
to the matrix Riccati equations (MREs) 



U' = Mq-UK^U, 
B' = -Nh+BK^B. 



(30) 
(31) 



Similar to the linearization of the scalar Riccati equations mentioned in J3J - I© °f the Introduction, the 
MREs l|5U )l . l|3T )l can be linearized by an Ansatz [231121] 



U 



VW~\ V,W € 



^2x2 



B = XV 
As result we arrive at the equation systems 



ITeC 2 



V 

w 







M 




V 

w 



X' 
V 



det(W) ^ 0, 
det(Y") ^ 0. 







Mi 




(32) 
(33) 



(34) 



The 4x2 matrices 



V 

w 



V 

w 



X 
Y 



VG 
WG 



are defined up to GL(2, 



X 
Y 



XG\ 
YGi 



x GL(2, C)— transformations 
G ,Gi e GA(2,C) 
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and can be interpreted as homogeneous coordinates of two points on a complex Grassmann manifold G2(C 4 ) 
which consists of 2— dimensional complex subspaces in C 4 (see, e.g. The matrices U = VW~ X and 

B = XY^ 1 are the corresponding affine coordinates of these points. 

Differentiating (|34[) and substituting V = KqW, X = —K\Y' it is easily seen that the equation systems 
(I34[l are equivalent to the second-order matrix differential equations 

(d r K d r - M ) W = 0, 

{d r K 1 d r -M 1 )Y = 0. (35) 

This implies that the matrices W — r -1 W ', Y = r~ x Y should be formal (non-normalized) solutions of the 
eigenvalue equations for the dynamo operator matrices Hi [ao], Hj-Jai], respectively 

H h [a ]W = EW, H h [ ai ]Y = EY. 

A comparison with the simple QM model from the Introduction shows that the intertwining operator matrix A 
should be expressible in terms of W or Y, and that W and Y should be connected by a product invariant like 
0. With the help of (JUJl, J22J and we find 

A = R(ip- Y'Y- 1 ) = (ip + KoW'W^Kg 1 ) R. 

In order to obtain the product invariant which connects W and Y, we use a slightly modified version of (|23|) 

U = RB^R- 1 -R(Rt)' 

and substitute from - J33I) 

U = RR i W'W-\ B* = -(Y i )- 1 (Y i yR*R 

so that 

W'W- 1 = -(Rt)- 1 (Y i )- 1 (Y i )'R i - {R 9 )- 1 ^)'. (36) 

This equation is of the type g — gin, {d r g)g~ x = (d r gi)gi + gi(d r n)n~ 1 g^ 1 . Hence, integration of 1|36|) yields 
the product invariant 

Y^R^W = C, det(C) ^ 0, 

with C a constant non-singular matrix. 

So far, we have obtained a 1 : 1 generalization of the intertwining technique from the simple QM example 
described in the Introduction to our J— symmetric dynamo operator model. It remains to test whether the 
MREs of this model are consistent. This will be the subject of the next section. 

4 No-go theorem 

In order to test the pair of MREs l|3U|) , l|31|l for consistency, we make use of (|23|l , (12411 as well as the relation 

N + K^N^K-l = K^K[ = K[ 

and transform the MRE for U [eq. I|30|l ] into an equivalent MRE for B. As result, we arrive at the following 
pair of MREs 

B' = R^MoR - K^BB + BK[ + [NN* + (iV")'] Ki, (37) 
B' = -Ah + BK^B (38) 

which should be satisfied simultaneously. The corresponding consistency test will be performed in two steps: 

1. From the limiting behavior at r — > we will derive a relation between a- n( i h- 

2. We will extract from eqs. iJSTJ, l|38() a system of non-linear ODEs for the helical turbulence functions 
ao, Oi\ and for the components b±, . . . , 64 of the matrix B. By mutual substitutions of these ODEs we will 
find an inconsistency which can be interpreted as a no-go theorem. 



|3Jl Limiting behavior at r — > 

From the assumed non-singular behavior of the helical turbulence functions at r — ► follows that they can be 
approximated as 

ao,i( r -»• 0) w c ,i + 0q,i7" + C(r 2 ) , c ,i ^ 0. 
Substituting this approximation in a slightly rewritten version of the defining equation l|35() for the matrix Y 



h(h + l) T ( E 



-ai 



oti E — a( 



Y = (39) 



Z^y'y- 1 m -iir- x /+ — a_ +C(r) , (40) 



we obtain the estimate 

K(r ->())» r-' 1 (/ + y CT _r + G(r 2 )) (r 2/l+1 C+ + C_) , 
where C+, C_ are arbitrary non-singular constant matrices det(C±) =/= 0. Correspondingly, it holds 

S = -X 1 y'y- 1 « Zxr- 1 /- [c 1 / 1 r- 1 + (? 1 + l/2)ai]a_+C(r) . (41) 
Comparison of (|41() with i|22|) shows that the components 62 and 64 of the matrix B vanish at least as 

b 2 ,b A ~0{r) forr^O. 

Furthermore, we find with the help of eqs. I|25l) . I|2t)[l and (|27|) that q ps 0{r) and, hence, ao/co = ai/ci, as well 
as <?', /, /' « O(l) what implies N, N$, (JV&)' « 0(1). 

We are now well prepared to perform a partial consistency test of (|37|) and (|38(l by comparing the singular 
terms of these equations in the vicinity of the origin r = 0. From the MREs (|37|l and (|38() we find 

-^K'iZ-Z' = + x ) j _ K^ZKiZ - ZK[ + 0(1), (42) 

-Kl x K' x Z-Z' = - h{h r ^ 1) ' l + ZZ + 0(l), (43) 
respectively. Substituting Z from (|40|l and equating the coefficients of the r~ 2 , r _1 — terms we obtain from eq. 

Zi = Zo + lj 01 = 

and, hence, also ao = 0. Eq. I|43|l is automatically satisfied, because Y is defined by the corresponding linearized 
equation l|39l) . The incremental relation l\ = lo + 1 is well known from ladder operator constructions for 
spherically symmetric Hamiltonians in QM This is not surprising, because this ladder operator construction 
can be recovered from the intertwining construction (|10() for the a 2 — dynamo operator matrices by the two-step 
transition: 1. ao = ol\ = a, 2. a — > 0. 

|4j 2 Systems of coupled non-linear ODEs and their inconsistency 

The system of eight coupled non- linear ODEs for the components b\ , . . . , 64 of the matrix B is easily obtained 
from the MREs J§ZJ|, Ipffil. e.g. with the help of the matrix multiplication package of M ATHEM ATICA© . For our 
analysis it is sufficient to consider only the simplest four equations of this system, i.e. the cr + and / projections 
of gJJ and 

b' 2 = 26i6 2 + ai(l + 6 2 ) (44) 

2 

= -2M2-— (45) 

b[ = b\ + b 2 b 3 — b\ — E— h ^ x 9 + ^ + aibib 2 (46) 

= -b\ b 2 b 3 + bl + E+ lo(l ° + 1} - a[b 2 + ( 3- + q '- q \ (47) 

2 

Equating the right-hand-sides of (|44|l . (|45|l and using 62 = 2g/ai from l|27|) we are able to express 61 as 

&1 = - 4g2+ Q Q ° + Qx . (48) 
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Taking into account that q = d r ln(ao/ai)/2 according to and that the helical turbulence functions «o and 
ct\ do not depend on Iq or l\ we conclude from equation (|48|) that 61 should not depend on Iq or l\ too. On the 
other hand, addition of (|46|l and (|47|l together with the relation Iq = li — 1 gives 



, 21 1 / ex. 1 \ ' t j 1 

26' 1 = -^+2g 6 X -^ 



<? 2 



what by integration leads to a function 61 which depends on 1%. I.e. the term depending on l\ cannot be 
compensated by a combination of ^—independent terms. This is an obvious contradiction to (|4H|) and we have 
to conclude that the consistency conditions - l|18|) cannot be fulfilled simultaneously This means that we 
are lead to the 
No-go theorem: 

The structure of the MHD a 2 — dynamo operator matrix is incompatible with an operator intertwining technique 
which is based on first-order differential intertwining operators. 

A similar situation occurs also for three-dimensional spherically symmetric models in QM J5|. There the 
I— dependent centrifugal term sets so strong restrictions on the form of the allowed potential that an intertwining 
construction built on first-order differential intertwining operators is only possible for the following three cases: 
the constant potential V(r) = const, the Coulomb potential V(r) oc 1/r, and the potential of the three- 
dimensional isotropic harmonic oscillator with V(r) oc r 2 . Richer classes of allowed potentials are only found 
for models in their s states, whe n I = 0. Such states are a priori excluded for the a 2 — dynamo operator matrix 
due to its construction [see Eq. l)A5jl ]. 



5 Concluding remarks 

In the present paper, we have tested the MHD a 2 — dynamo operator matrix for its compatibility with the 
simplest variant of an intertwining construction based on first-order differential intertwining operators. The 
operators have been chosen in accordance with the fundamental J— symmetry (pseudo-Hermiticity) of the 
operator matrix and lead to a set of six matrix equations as consistency conditions. With the help of a step- 
by-step reduction of the complexity we have extracted their basic structural elements and have shown that 
they contain an intrinsic inconsistency. So, we have to conclude that the structure of the a 2 — dynamo operator 
matrix is not compatible with the considered first-order differential intertwining Ansatz. This fact is the subject 
of the formulated no-go theorem. 

It remains to test whether intertwining constructions can be built from second-order or higher-order differen- 
tial intertwining operators. Energy shift operators based on second-order differential expressions are known for 
harmonic oscillators with time-dependent frequencies and additional 1/r 2 — term |25| as well as for the spherically 
symmetric oscillator and the Coulomb potential A generalization of the technique to the MHD a 2 — dynamo 
operator matrix seems realistic. 

Another approach for a clarification of the considered isospectrality problem could consist in a generaliza- 
tion of the Gelfand-Levitan technique for vector- valued Sturm-Liouville problems [2fi| . Concerning its general 

structure, the a 2 — dynamo operator matrix -ff/[a] is a singular non-self-adjoint matrix Sturm-Liouville operator 
which by a unitary transformation can be recast into the standard form 

-d r P 2 {r)d r + P (r) . 

In 1998, Jodeit and Levitan |2(i| analyzed the isospectrality problem for matrix Sturm-Liouville operators with 
P2M = I and Po(r) a symmetric matrix. They showed that if two vector- valued Sturm-Liouville problems 
are isospectral then the eigenfunctions of one problem can be constructed from the eigenfunctions of the other 
problem with the help of a matrix Gelfand-Levitan transformation. So, a generalization of this technique to 
Sturm-Liouville problems with non-symmetric Pq(t) and -P2M 7^ I2 would naturally cover the isospectrality 
problem for the MHD a 2 — dynamo operator matrix. 
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A Derivation of the a 2 — dynamo operator matrix from the mean- 
field induction equation 

For completeness we sketch here the main steps of the derivation of the 2x2 operator matrix Hi [a] for a model 
with helical turbulence function a(r). The outline follows the technique for models with a = const as presented 
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The spherical MHD mean-held a —dynamo in its kinematic regime is described by the induction equation 
for the magnetic held 

ftB = Vx (oiB) + z/ ro AB (Al) 

supplemented by the condition V-B = 0. The magnetic diffusivity v m is assumed to be constant and the helical 
turbulence function a to depend only on the distance from the origin a — a(r). Decomposition into toroidal 
and poloidal components B = B t + B p and setting B p = VxAt allows for a decomposition of the induction 
equation IjAljl 



d t B t = Vx(a Vx A t ) - u m Vx VxB t 
dt A t = aB t - v m V X V X A t . 



(A2) 
(A3) 



Furthermore, the fields Bt and At can be represented as 

A t = -r x V Fx, B t = -r x VF 2 , 
where F± and F2 are single- valued scalar functions which are normalized on the unit sphere S 2 by the condition 



Pi 2 du = 0. 



(A4) 



r x V APi 
'1 



r XV 



(9 r a)(5 r rPi) + aAPi 



With the help of the relations 

A(r XVPi) 
Vx [a Vx (-r X VP 1 !)] 
ar X VP 2 

equations (|A2|) and HA3|> can be rewritten as 

r XV [u m AF 1 + aF 2 -d t F 1 



r x V(aF 2 ) 



r X V 



f m Af 2 - i(a r a)(9 r rP 1 ) - aAF 1 - d t F 2 



= 0, 
= 0. 



It follows that the expressions in the square brackets are functions of r and t alone which must vanish due to 
the normalization condition 1A4|I and its implication J g2 AFi^du = 0. By re-scaling of r and t one sets the 
magnetic diffusivity to unity v m — 1 and the boundary conditions at r = 1. 
With the help of a series expansion in spherical harmonics 

F 1)2 = ^ e^» t ^f n) (r)y ( ra (M) £ i 2 (fl,r 2 *) ® L 2 (S 2 ,^), O = [0, 1] 

l.m.n 



one obtains the eigenvalue problem 



AiF (l,m,n) + a p(l,m,n) 



— ^i, n F 1 



(Z,m,n) 



A; P. 



(j,m,n) _ i(gi ra )(a r rPf ' ro ' n) ) - a&iF[ l > m > n) 



,ri ± 2 



Here we used the notation A; = \d r r 2 d r — and the fact that due to the symmetry of the dynamo 

configuration the eigenvalues \i :U depend only on / and n. We note that the normalization condition l)A4|l 
implies 

(A5) 



p(l=Q,m,n) _ q 



Finally, the substitutions p = —i(d r + 1/r), -01,2 = F : 



(Z,m,n) 
1.2 



€ L 2 (il,r 2 dr) lead to the eigenvalue problem for 



the a 2 — dynamo operator matrix P/;[a] as it is given in eq. (JHJ of the Introduction. 
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